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We study the effects of Landau level mixing in the limit of weak electron interaction. We use a numerical 
method to obtain the two- and three-body corrections to quantum Hall pseudopotentials, which are exact to low- 
est order in the Landau level mixing parameter. Our results are in general agreement with certain analytic results 
(some derived here, some derived by other authors) in the thermodynamic limit. We find that the convergence to 
this thermodynamic limit can be slow. This suggests that errors could occur if one tries to use pseudopotentials 
derived in a thermodynamic limit for numerical work on finite systems. 



I. INTRODUCTION 

Fractional quantum Hall effect is typically studied in the 
limit where the cyclotron energy hu> c is large, so one can fo- 
cus attention on only a single Landau level (LL). However, in 
typical experiments, the interaction energy scale -E C ou/oro6 = 
e 2 /(e£o), with Iq the magnetic length, is often on the same 
order as the cyclotron energy, thus making Landau level mix- 
ing an important consideration. Historically, the difficulty of 
including LL mixing into calculations, along with the general 
belief that LL mixing would change predictions quantitatively, 
but not qualitatively, has generally prevented researchers from 
addressing this issue in much detail. Relatively recently, how- 
ever, it became clear that one could not always ignore LL mix- 
ing so easily. Crucially, it was pointed out that LL mixing 
would be required in order to break the symmetry between 
two possible (nonabelian) quantum Hall states that might be 
realized at v = 5/2 — the Moore-Read^ state, and the so- 
called AntiPfaffian 2 . In this context, the current authors used 
a Landau-level truncation scheme to conclude that in typical 
samples the Anti-Pfaffian is favored^. 

Although the experiments are not usually in the limit of 
small LL mixing, it is nonetheless natural to consider the case 
where LL mixing is weak. To this end, Bishara and Nayak 4 
used a systematic perturbative approach to obtain corrections 
to leading order in the parameter k — E cou i om b/hhj c . As 
pointed out in that work, the effect at leading order can be 
summarized as being a correction to the two-body 5 and three- 
body^ pseudopotentials. These corrections were used in an 
exact diagonalization study in Ref.[U 

Almost two decades earlier, a seminumerical LL mixing 
study was undertaken by Rezayi and HaldaneP to obtain the 
leading perturbative corrections to the two-body interaction. 
Aware that the three-body interactions would be crucial, in 
the process of writing Ref. 3, the current authors explored ex- 
tending the numerical technique of Ref.|9]to obtain three-body 
corrections as well — the results of which will be presented 
in the current paper (these results also update the results of 
Ref. |9] which were necessarily rough given the limited com- 
putational ability of the time). Unfortunately, we found that at 
least some of our results disagreed with the published results 
of Bishara and Nayak 4 . We later supported our numerical re- 
sults by an analytical calculation valid in the lowest LL, which 
we present in the appendix of this paper. 

In discussions with Nayak as well as with Matthias Troyer, 



it was suggested that the disagreement in results was due to 
a normal-ordering error in Ref. 4 (This is now discussed in 
some detail in Ref.fTOb. However, even fixing this error it was 
not obvious that the analytical results entirely agreed with the 
numerical calculation. 

More recently another fully analytical calculation was per- 
formed by Sodemann and MacdonalcPil While the results 
agreed with our numerical calculation in the lowest LL (LLL), 
some disagreements remained in the first exicted LL (LL1). 

Still more recently, Peterson and Nayak have produced a 
new manuscript 10 performing an analytical calculation which 
now agrees with our numerical results (their manuscript also 
discusses finite thickness corrections and graphene, as well 
as the conventional pure Coulomb interaction). Since these 
results are controversial, and since it is obviously very easy 
to obtain incorrect results, we have decided to publish our 
numerical results to support the conclusions of Peterson and 
Nayak. In so doing, we will also point out a key potential lim- 
itation of the analytic techniques of these references — that 
they apply only in the thermodynamic limit. 



II. METHOD 

The method roughly follows that of Ref. 9 We work on a fi- 
nite sphere with a monopole of strength N<j, flux. We start with 
some number of completely filled Landau levels. We would 
like to consider the interactions of n = 2 or 3 additional elec- 
trons added to the next empty LL, which we call the valence 
LL. For the two (three) particle case we consider a basis of 
states with two (three) electrons added to the valence LL. For 
simplicity we consider only states with total L z = 0, which 
will be sufficient for our purposes. One can diagonalize the 
L 2 operator within this space to pick out states with angular 
momentum L which we call \L, n) where n = 2, 3 indicates 
whether we are considering two or three electrons added to 
the valence LL. While specifying L picks out a unique wave- 
function for the case of n — 2 electrons added to the valence 
LL, for the case of n — 3 electrons added to the valence LL, 
the quantum number L alone does not uniquely specify the 
state for L = 9 and L > 11, in the spin polarized case, and 
for L > 4 in the unpolarized case^El. In these cases, we need 
additional quantum numbers to specify a specific wavefunc- 
tion. Although we will not delve into such complexity in the 
current paper, the described technique can easily be extended 
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to handle these cases as well. 

Our Hamiltonian can be abbreviated as 

H = K + U 

with K the kinetic terms and U the interaction term which we 
assume to be purely two-body (and below we will take this to 
be of the Coulomb form). The above mentioned \L, n) states 
all have the same kinetic energy K® 

K° n \L,n) = K\L,n) 
We next apply the interaction to our basis of states 

\L\n) = U\L,n) 

It should be noted that these wavefunctions contain compo- 
nents with electrons above the valence LL, as well as holes 
below the valence LL. Crucially, for numerical work we must 
truncate the Hilbert space that we will consider by not al- 
lowing electrons to occupy Landau levels above some level 
LL max . Note that written in the standard single electron or- 
bital basis, the number of components contained in a one of 

these wavefunctions \L, n) is only polynomially large in 
and in LL max . 

Note that the Coulomb energy of these wavefunctions with- 
out considering LL mixing can be written as 

V£ >n = (L,n\U\L,n) = (L,n\L^). 

For the case of n = 2 electrons in a valence LL, these val- 
ues are precisely the pseudopotentials in the absence of LL 
mixing. They are proportional to E cou i om b. 

The leading LL mixing correction to the energy of the \L) 
state is then given by 

SV L , n = ^\P—?—,P\L^) 

where P projects to the space where K 7^ K®. In the case 
of n = 2 electron SVl,2 is precisely the leading LL mixing 
correction to the V® 2 pseudopotentials. These energies are 

proportional to nE CO uiomb = (e 2 /e£o) 2 /Swc- Note that in the 
Results section below we will always quote energies in units 

of kJL coulomb- 

For the n = 3 electron case, we want to find the "irre- 
ducible" part of the three-body interaction — that is, the part 
of the energy shift that cannot be accounted for by a correc- 
tion to the two-body pseudopotential. We first calcualte V® 2 
and SVl. 2- We then consider a modified Coulomb interaction 
within the single valence LL that consists only of V® 2 + SVl,2 
and contains no LL mixing terms. We then evaluate the shift 
in V® 3 resulting from having added SVl,2 to the t> are pseu- 
dopotential V® 2 - This shift is subtracted from SVl 3 to give 
the irreducible three body pseudopotential. Note that obtain- 
ing the three-body term thus requires obtaining the two-body 
term first. 

The described technique determines the leading order LL- 
mixing corrections to the interaction. However, we must recall 
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Table I. Predictions for the two-body pseudopotential leading correc- 
tions due to LL mixing for electrons in a spin polarized LL. The top 
table is results for the LLL, the bottom table is the results for LL1. 
Again, we quote differences in pseudopotentials to remove physi- 
cally unimportant self-energy terms. As suggested by the bottom 
panel of Fig. [T] in the case of LL1, one should assume extrapolation 
error of possibly as much as .005. As in that figure the extrapolation 
considers LLmax < 6 and < 20. 



that we have truncated the Hilbert space to a finite number of 
LL's and we must numerically extrapolate to infinite LL max . 
Finally, since we are interested in comparing to results de- 
rived for large (planar) systems, we must also extrapolate to 
the limit of large corresponding to a large sphere. Thus, 
the numerical results presented below are the results of a dou- 
ble extrapolation. Note however if pseudopotentials are to be 
used on finite sized systems, for consistency, it is preferable 
to use pseudopotentials explicitly derived for the system size 
in question. We will discuss this issue further below. 

HI. RESULTS 

All numerical results will be quoted in units of nE cou iomb- 
To avoid having to worry about electron self-energies (which 
may converge more slowly) it is more convenient to keep track 
of differences in pseudopotential coefficients. This is natu- 
ral given that adding a common constant to all pseudopoten- 
tial coefficients gives an overall shift in energy but leaves the 
eigenstates of the Hamiltonian unchanged. To calibrate our 
technique, we consider the two body pseudopotential in the 
lowest LL (LLL) as shown in the top of Fig. [TJ Analytic re- 
sults are derived for this case in the appendix (see also Refs.fTOl 
and ITTb . In this case, the numerical extrapolation is fairly 
clean. However, in the case of the two-body corrections in 
the first excited LL (LL1), as shown in the bottom of Fig. [TJ 
the extrapolation involves a guess at the proper extrapolation 
functional form; as a result one should assume a bit more er- 
ror in the result. Nonetheless, comparing our result with the 
analytic prediction of Ref. 10 still obtains agreement within 
about 1%. 

In Table [I] we show LL mixing corrections to the two-body 
pseudopotentials in the LLL and in LL1 . In the LLL the agree- 
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Figure 1. (color online) Sample Double Extrapolations. The multiple 
curves in the main plots correspond to different N#. The intercepts 
(infinite LLmax) are then plotted in the inset and extrapolated to in- 
finite N^. Top: Extrapolation of the two body LL-mixing correction 
terms 8Vz — SVi for electrons in the LLL. Here LL max ranges from 
2 to 6. iV^ ranges from 12 to 24. The double extrapolation yields 
a value 0.0273. compared to the analytically exact result 0.0272267 
(See Table I, Refs 10 and 1 1 and the appendix below in this paper). 
Bottom: Extrapolation of the two body LL-mixing correction terms 
SV$ — SVi for electrons in the 1st excited LL. Here LL m ax ranges 
from 3 to 6. N$ ranges from 12 to 20. Due to the curvature in these 
plots, one must assume a bit more error on the extrapolation. The 
double extrapolation yields a value 0.1089. The analytical reslult of 
reference 10 yields .1104. This seems well within the range of pos- 
sible extrapolation error. 



LLL 


Numerical 


Analytic^ 


V B - V 3 


.0214 


.0214 


V 6 - Vi 


.0060 


.0074 


v 7 - v 3 


.0239 


.0240 


V 8 - Vi 


.0123 


.0133 


LL1 


Numerical 


Analytic™ 


V s - Vi 


.0109 


.0093 


V 6 -Vi 


.0027 


.0048 


V7-V3 


.0191 


.0152 


V a -Vi 


.0177 


.0138 



Table II. Predictions for the three-body pseudopotential coefficients 
due to LL mixing for electrons in a spin polarized LL. Analytic re- 
sults are taken from both Refs. [To] and [TT] Again, we quote dif- 
ferences in pseudopotentials to remove physically unimportant self- 
energy terms. The numerical results in the LLL (top) are extrapolate 
extremely accurately to the thermodynamic limit. However, for LL1 
(bottom) the agreement is less good. (See text for discussion). 



ment between the numerics and the analtyics is excellent. In 
LL1, the agreement is less precise, but considering the diffi- 
culty in extrapolation (See Fig. [T}, the agreement can still be 
considered to be fairly good. 

In Table|H|we show the three-body pseudopotentials result- 
ing from LL mixing for a spin polarized valence LL. In the 
LLL the agreement between numerics (this work) and analyt- 
ics (Refs.fTOlandfTlTl is quite good, whereas in LL1, the agree- 
ment is much less good. To examine this situation a bit more 
closely, we show an example of some of the extrapolations in 
Fig. [2] Note that there is quite a bit of change in the value of 
the pseudopotential as a function of N^, which makes extrap- 
olation difficult. Nonetheless, given the uncertainties in the 
extrapolation, the agreement with the analytic results appears 
acceptable. 

The fact that the convergence with system size is somewhat 
worse in LL1 than in LLL is to be expected being that all 
correlation lengths are larger in LL1 (the single particle wave- 
functions are in some sense "fatter"). In addition, working 
in LL1 is far more computationally demanding than in LLL 
since one needs to account for a far greater variety of exci- 
tations that involve exciting electrons out of the fully filled 
Landau levels (including both spin species). This has limited 
our numerical calculations in LL1 to somewhat smaller sys- 
tems. While these particular calculations were performed on 
a not-very-powerful computer, the results are clear enough to 
demonstrate the general trends, and also to confirm agreement 
with analytic work by other authors. 

The slow convergence of the three-body potential raises a 
serious issue. To a large extent, the purpose of deriving pseu- 
dopotentials is to use them for numerical exact diagonaliza- 
tions for finite systems. However, there is a clear inconsis- 
tency in using pseudopotentials derived in the thermodynamic 
limit for exact diagonalizations of a finite sized system. While 
the size dependence of the two-body pseudopotentials is rela- 
tively minor, the size dependence of the three-body term ap- 
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Figure 2. (Color online) Sample double extrapolations of the three 
body pseudopotential Vs — Vz (top) and Ve — Vz (bottom) in LL1. 
The extrapolated values are 0.109 and .0027 respectively. Note that 
there is quite a bit of slope with respect to N^, which makes accurate 
extrapolation difficult. Secondly note that the extrapolated result is 
close to zero on the scale of the values measured at finite N$ which 
also makes the extrapolation difficult. 



pears to be very significant. This is also true for the LLL as 
well (perhaps to a slightly lesser extent), as shown in Fig. [3] 
This raises the possibility that, for cases where the three-body 
potential is important in determining the physics, one will not 
see the true thermodynamic behavior until very large system 
size. We suggest that finite size extrapolations may have fewer 
errors if diagonalizations of a given size system use pseudopo- 
tentials designed for that particular size. 

To elaborate on the issue of extrapolation, we recall the 
work of Ref. [12] There it was pointed out that extrapolation 
converges rapidly if (a) one uses interaction (such as chord 
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Figure 3. (Color online) Sample double extrapolations of the three 
body pseudopotential V5 — V3 in the LLL. The extrapolated values is 
.0214 Note that there the value of the pseudopotential taken at finite 
Nff, can be quite different from the thermodynamic limit. 



distance) appropriate for each particular size system and (b) 
one rescales the interaction strength to keep the charge den- 
sity constant. For this reason, in the history of quantum Hall 
numerics, the vast majority of calculations have used pseu- 
dopotentials derived for each particular system size in ques- 
tion. While we cannot say with certainty that other schemes 
may not be made to work, it seems very risky to throw out the 
inherent consistency of this established approach. 



IV. SUMMARY 

The LL mixing corrections presented here should be of 
value in the analysis of physical systems where LL mixing 
terms are weak (very high density quantum Hall samples). 
Moreover, the results here have been used to independently 
establish the validity of analytical results of Ref.|T0]and some 
of the results of Ref. 1 1 Finally, we point out a potential in- 
consistency in using pseudopotentials derived in the thermo- 
dynamic limit for calculations on finite systems. 

After the initial submission of this paper, a revised version 
of Ref. 1 1 appeared which is in full agreement with this paper 
as well as with the results of Ref.flOl 
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Appendix In this appendix we analytically calculate the 
leading LL mixing correction to the two-body pseudopoten- 
tial in the lowest LL. 

For a single electron in a magnetic field defineP 



a = V2 (lB z + 

b = V2 (id z + ^z 



Where £ is the magnetic length. Noting that d\ 
z' — z we have 



[o,ot] = [6,6t] = i 



The Hamiltonian is 



H = hu c (a) a + — ) 

and the z-component of angular momentum is 

L z = h(tfb-Or~a) 



(1) 

(2) 
-d z and 

(3) 
(4) 
(5) 



The single particle eigenstates with energy (n + \)?vjj c and 
L z = H(m — n) can be generated as 



Km) - SS^TO 

v n'.mi 



(6) 



where |0, 0) is the fiducial state for which a|00) = 6|00) = 0, 
which takes the form 



(rl0,0) 



1 



-\z\ 2 /4( 2 



(7) 



The explicit form of all the single particle eigenstates can be 
shown to be 



(r\n, m 
(-1)" 



(8) 



\zl 
2£ 2 



-\z\ 2 /U 2 



(9) 



Now consider two electrons with coordinates z\ and z^. We 
can transform to center of mass zq and relative coordinates zr 



zc 



1 -1 

1 l 

2 2 



Z\ 
Zl 



(10) 



This transformation has unit Jacobian. Correspondingly we 
have 



d z 



/ I _I 

_ | 2 2 
1 1 



(11) 



which maintains the desired commutations [d ZR , zr] — 
[d zc ,z c ] = 1. 

It is convenient to write out the reverse transformations as 
well 



Zl 

zi 



and 



3 1 
-3 1 



-1 h 



zr 
zc 



Substituting this into our expression for <2j we obtain 



Ol,2 = V2 



1 



1 



(12) 



(13) 



(14) 



Defining £r = £\[2 and £c = £/\2 we can rewrite this as 



Ol,2 



V2 



(£ c d zn + -i 



1 

Ur 



zc 



± V2 ( £ R B ZR + —z R 



= -7^(ac±a R ) 



where we have defined the operators 



a R = V2 \ £ R d ZR + -^ZR 
a c = V2 \ £cd Zc + 



(15) 
(16) 

(17) 
(18) 



which correctly satisfy [a R , a R ] — [a c , (r c \ = 1. 
Similarly we may define operators 



''/; = v 2 ( £ R d ZR + ^-z R 



b c = V2[ £ C d. 



1 

Mr 



-zc 



such that 



& i,2 = ~^( b c ± h R) 



(19) 
(20) 

(21) 



and again [b R , b R ] = [be, bg]=l. 

In terms of these new variables, we have the Hamiltonian 



H = huj c ya\ai + 0^0,2 J (22) 
= Huj c (a R a R + a c a c ^J (23) 
And similarly the total z-component of angular momentum 



L z = h 



(b R b R +bl 



c b c - a' R a R - a c a c 



(24) 
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The key here is that a radial interaction V(\ri — f2|) = 
V(|z.r|) between the two particles cannot effect the center of 
mass degree of freedom, so we may freely choose all of the 
ac and be oscillators to remain in their ground states. Our 
kets will then be of the form \N, M) where N and M are the 
occupancies of the a R and b R oscillators. 

Our initial state will have both electrons in the LLL and will 
give them relative angular momentum M. Thus we choose an 



initial state |0, M) for the a R and b R oscillators. A rotation- 
ally symmetric interaction cannot change L z — b R b R — a R a R , 
thus we can only have matrix elements of the form 

(N,M + N\V\0,M) (25) 

We can directly evaluate this matrix element as 



11 (J*.\ M rMW^-WV^ 



1 2tt£ 2 r (M)\ 1 


( ZR \ 







(26) 



AH 



2 \ M 

rdrV(r)l") L ff(^)e^^ (27) 



Let us now plug in the coulomb form of the interaction V = e / (er) and making the substitution x = r /2£ R we obtain 



(N,M + N\V\Q,M) = ^\J—^-— I dxx M - 1 / 2 L^(x)e- x (28) 



e£ R V2 \ Ml(M + N)\ J 



The integral is standard on integral tables and can be per- 
formed to give 



{N,M + N\V\0,M) 



e 2 T(M+ 1/2)T(A + 1/2) 
2yJnN\ M\ (A + M)l 



(29) 



and note that this is now written in terms of the original mag- 
netic length. 

As a check we can determine the Coulomb pseudopoten- 
tials in the LLL with no LLL mixing. These would be given 
by 



V M = (Q,M\V\0,M) 



e 2 T(M + l/2) 
H MY2 



(30) 



which matches the well known expressiorP^. 

Now in second order perturbation theory the change in en- 
ergy of the |0, M) state should be given by 



5V M - J2 



\(N,M + N\V\0,M}\< 



N>0 



-huj c N 



(31) 



This result agrees with the analytic expression given in 
Refs.QlandQT] 
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